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ABSTRACT 


The phase-plane method is a graphical method for 
linear and non-linear second-order systems. There are 
many techniques for obtaining the phase portrait which 
consists of a number of phase trajectories on the phase- 
plane. From the summary and discussion the isocline 
method is a most general and useful method. The only 
difficulty is in the labor required. Solution by digi- 
tal computer reduces this difficulty and makes the iso- 
cline method a much more useful method in non-linear 


systems analysis and design applications. 
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I. INTRODUCTION 


JEN CONDITIONS FOR LINEARITY AND DEFINITION 


OF NON-LINEARITY 


It is useful to define linearity and then any systen, 
component, or phenomenon, which does not satisfy this de- 
finition will be considered non-linear. The simple mean- 
ing of linearity is that if the relationship between output 
y and input x is plotted the result is a straight line. Or 
that if the output is proportional to the input the process 
must be linear. So we can simply say the essential condi- 
tion for linearity is that the law of superposition applies, 
and linearity is a specification which therefore limits the 
field of activity. Non-linearity is a "non-specification" 
and the field is unbounded, linearity being a particular 


feet of it. 


oui CLASSES OF NON-LINEARITY 

Non-linearity in physical systems cannot be avoided, 
and in fact can sometimes be used to improve the system's 
performance. The non-linearities encountered in control 
systems may be classified as accidental non-linearities 
and intentional non-linearities. The accidental class con- 
sists of those non-linearities which are inherent in the 
components used in the system. Some of these are coulomb 
friction, stiction, backlash, binding, dead zone, satura- 


tion, etc. Their existence is accidental in the sense that 


15 


they are not deliberately inserted in the system. Inten- 
tional non-linearities are those which are introduced de- 
liberately, either to obtain desired performance features 
or for weight, space, or economic advantages. 

Clearly the simplest class of non-linearity is the in- 
stantaneous, a member of which only changes its input u(t) 
to some output x(t) without any delay, the ratio x(t) /u(t) 
being.a function of u(t). Thus x(t) = flu(t) ] and it can 
be represented in schematic form (see Fig.(1-1)), having a 
transfer function f(z)/z which is independent of time. 
Saturating amplifiers and relays are examples of instan- 
taneous non-linearities. In practice, of course, such a 
non-linearity will almost always be associated with some 
linear (or non-linear) weighting function resulting from 


inertia and damping as in Fig. (1-2),so that 
(SP Seu) 4 (1-1) 


In general, the instantaneous non-linear characteristic is 
combined with inertia and damping in the same component, in 


which case it cannot be separated as in Fig. (1-2). 


x(t) 





E2) which 


is independent of time. 


Fig. l-l. A transfer function 
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Fig. 1-2. The instantaneous non-linear 


characteristic is combined with 
inertia and damping cannot be 
separated as in figure. 


By convolution integral 


c 


ole = | Wihit-tidr , and 
O 


changing the variable gives 
e 


y(t) = [| x(t-7,)h(7,)ar, 
O 


Assume 


V7 


Then 


oc 


mre) | x(t-T,)h(7,) dr, 
O 


| f[u(t-7,) Jh(r,)dr, 
O 


II 


E 
I K[u(t-7,),7, lary (1-2) 
- © 


K[u(t-7,),7,] is called the Kernel function and is seen to 
depend on both the dynamics of the system, £(z)/z and h(t), 
and the characteristics of the input process, u(t). 

Zadeh has classified non-linear systems on the basis 
of multiway integrals and Kernel functions. Thus those 
that can be represented by a single-way integral belong to 
the first class, Ny which therefore includes the subclass 
of all linear systems and also the degenerate zero-order 
class of instantaneous non-linearities. 

Then the first-order non-linear differential equation 
is 


= + F(x) = F, u(t) J 


Os 


ax 
spe 1 eB) J 8) 
Integrating gives 


c 


ee | K{u(t-7),7T]d7 
O 
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Hence, the second-order non-linear differential equation is 


£, (244 fi (  )4 ST dans geile 


f(x, x, x) = uid. 


Integrating twice gives 


mia ae | ( K[u(t-7,), u(t-7,),7,,7,]4T,aT,. 
O O 


So, the second class, Nos contains all those for which 
the input/output relationship can be expressed as a 2-way 
integral and this includes the degenerate cases belonging 
to class n° 

The nth class, nA? contains all those systems for 
which the input/output relationship is a n-way integral, 
which means that they can be represented by nth order non- 
linear differential equations. 

Zadeh's classification is an important property in 


considering random processes in the presence of noise. 


oF NON-LINEAR ANALYSIS 

Many methods used in the analysis of non-linear prob- 
lems are effective only for limited types of problems. A 
few of them have been developed which have more general 


applicability to non-linear control systems: 


1S) 


1. Linear approximations. 

2. Graphical method - The phase-plane. 
3. Piecewise linear approximation. 

4, Numerical methods. 

5. *DESeri bingssunc cron, 


6. The analog/digital computer. 


D. PHASE-PLANE ANALYSIS OF NONLINEAR SYSTEMS 

The phase-plane method is a graphical method. There 
are many techniques for obtaining the phase portrait which 
consists of a number of phase trajectories on the x vs x 
plane or in some plane related to this plane through a con- 
venient transformation of variables. After the phase por- 
trait is obtained, by means of it, analysis of the follow- 
ing aspects and items is possible: 

1. Computation of the exact transient response for a 
given set of initial conditions. 

2. Prediction of the type of transient response 
(oscillatory or monotonic) with a certain range of initial 
conditions. 

3. Estimates of the amount of overshoot in response 
to a step input of given magnitude. 

4. Prediction of the system stability. 

5. Prediction of the existence and/or amplitude of 
limit cycles. 

6. Prediction of the existence of various modes of 


operation. 
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7. Study of the effect of a given type of non- 
linearity on the performance of the system as compared 
with a linear system, with a system having a different 
type of non-linearity, or with a system which is non-linear 
in the same sense but in different degree. 

8. Attempts to answer questions such as: 

a. Does this non-linearity improve performance? 

b. What changes can be made to suppress the limit 
cycle? 

c. How can we adjust this system to make it meet 
specifications? 

As stated above, we can understand a graphical method; 
it is usually relatively simple to utilize, and may be 
especially useful as an exploratory tool when a non-linear 
equation is first being attacked. There are many ways to 
obtain phase trajectories on the phase-plane. In the next 
chapter, the most useful methods will be summarized and 
compared. The advantages and disadvantages of the various 


methods will be pointed out. 


2. 


A. 


IT. SUMMARY OF USEFUL GRAPHICAL METHODS 
FOR CONSTRUCTING PHASE TRAJECTORIES 


ON THE PHASE-PLANE 


ISOCLINE METHOD 


1. Method-description 


The free motion of any second-order non-linear 


system can be described by an equation of the form: 


x + O(x,x)xX + V(x,x)x = 0 


where 
= glx 
x dt 
a 4 2 
site 
Let 
ax 
- — 2-1 
ae ( ) 
then 
——— &y =- @@,y)y — Vee (2-2) 


Equation (2-2) divided by (2-1) yields 


Let 


dy os (x + W(x x 
ax y . 
a = M where M is a constant. 


22 


Then 


My = @x,y)y + W(x,y)x (2-3) 


Equation (2-3) is the equation of isoclines. 

NOTES: a. If @(x,y) and U(x,y) are constants, then 
isoclines are straight lines. Otherwise, isoclines will 
be curves. 


b. For a first-order equation: 


Ox _ 


Isoclines will be £f(x,t) = M. 

c. Equations more than second-order need some 
transformation technique, so this method is good only for 
second-order systems. 

d. After obtaining the isoclines, phase trajec- 
tories are plotted on the phase-plane as follows: 

M,, Mo, M. --- represent the slope of each 
isocline. Start from same initial condition point, as in 
Fig. 2-1 point A. From point A draw dotted lines parallel 
to M 


M.- These dotted lines intersect the M, isocline 


ilk 
at a and be. Bis middle point of ab. Continuing in this 
wey, Obeain C, D ... points. Commect A, B 2.2 £0 Gbraun 


a phase trajectory from the given initial condition A on 


the phase-plane. 





2. Example 
Consider a viscous-damped servo 
ao. a0. 
Je + Fap t+ KO, = 0. 
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EF Ciiee 2 aL awe COnSELUGCELON “OLA 
trajectory using 


isoclines. 
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Dividing through by J: 




















2 
ante nS $. mA = 
ar J dt O 
SUbstTeuting 
a ; K 
a of 
and 
F 
c = <a 
2,/KI 
into equation (2-4): 
a°o a0. , 
+ 2cW zm +w"8 =0. 
ate cre n o 
Let 
dé 
ss 
at 
then 
2 
aa. a as ) ae ¢ ig 
dé. ao | dt Zz at 
sO 
waz 2 = 
Za 2CW Z + we Oo. = 0), 
Let 
GZ 
ap 
O 


Zo 


(2-4) 


Isocline equations will be 


- 2 
ZM = 2cu) Z Ww AO: (2-5) 


If c = 0.25 (linear servo system) 
and WwW = 1, 
n 
equation (2-5) will be ZM = -0.52z2 lc 
ti Mis equalyto 2,1,0, -0.5, =-l, -2 we@eecain seraigme 
radial lines (isoclines) as in Fig. 2-2. 


Assume an initial condition represented by A in”“the 


figure, namely an output position 8S = -2.7 and a position 
dé 
rate of change aa = 0, we can plot the phase trajectory as 


UE Ob Grille) Ar: 
3. Discussion 

a. From Fig. 2-2 it is seen that by plotting some 
isoclines on the Shasezpilene (more or less achenaiing on de- 
Sired accuracy) then any phase trajectory is easily plotted 
from given initial conditions on the phase-plane. 

b. By inspection of a few trajectories the sys- 
stem's behavior is determined. For example: In Fig. 2-2 
this system has damped° oscillations and is a stable system. 
Point B is an oscillatwen max. (6. =" | ery 

c. In this example, as applied to the case ofa 
linear servo system (C is constant), isoclines are straight 
lines, so it is simple. If C is not constant more labor 


1s required. 
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Fig. 2-2. Phase-plane plot of damped 
oscillation obtained by iso- 


cline method for a servomechanism. 


a 


B. | LIENARD METHOD 
1. Method Description 
Liénard's method is a means for drawing a direc- 
tion field in the phase-plane to guide the tracing of the 
trajectories of a differential equation. 

Using the isocline method, (see Fig. 2-3), draw the 
zero-Slope isocline on the phase-plane. Let P be any 
point of the plane at which the slope of the direction 
field is desired. ‘From P draw the line parallel to the 
Z-axis, which intersects the zero-slope isocline at Q. 

Let R be the projection of Q on the z-axis. Now, draw a 
short line segment through P, perpendicular to the straight 
line RP). Wier scamentersmtangene to the direction field 
of the integral curves of certain differential equations. 

In practice, graphical construction of the entire 
direction field is quite easy since the short segments at 
all the points on the extension of the line PQ can be drawn 
by using a compass Serle yo6 at she one point R. Note that 
this construction does, indeed, result in @® = O everywhere 
on the zero-slope isocline. 

In Fig. 2-3 the two angles @ and B are equal because 
of the perpendicularity of their sides at P and Q. The 
tangent of each angle can be evaluated in terms of the co- 
ordinates of the points P, Q and R and the ratio of lengths 
of unit distance along the two axes. This length ratio is 


defined by 


ARS: 












-Zero-slope 
isocline 


Fig. 2-3. Liénard's construction of 


phase-plane direction field. 
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_ unit distance for v 


7 unit distance stor Z (2G 
= cv 
aL ani FET ee ao (278) 


where M is the slope of the integral curves and/or of their 
direction field 
So) ae ees ¢ Se ©) 
The expression Zp = £(V,) is the relation defining the 
zero-Slope isocline. 


Since V. = V P is any point of the phase-plane. The 


Q P@ 
subscript can be dropped, giving 


I! 


tan p = = a (2-9) 


since 


tan w = tan B 


then 
Ean fo = KM = tan #8 = a ae 

m= 2 (2-10) 
K V 
dv 

" dz 

dz av saz az 
Mdt 7 dz dt 7 5,2 ls 
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Combining equation (2-10) and (2-11) 





az 7 
K wae ae ) We 6 (ey 


Any servo which obeys a differential equation of this 
form has a phase-plane representation with a direction 
field which can be constructed by Liénard's method. 

Note that the coefficient of Z2 is unity and that K 
specifies the ratio of vertical-to-horizontal unit distan- 


ces on the phase-plane. The zero-slope isocline is found 





2 
directly from the differential equation by setting e 5 
Gite 
equal to zero, since the isocline is the line for which 
| 2 
M= 0. So equation (2-11) gives o% = OF 
dt 


per lapses ty OWES 


Using the same example as in A, a linear servo- 
mechanism having a constant damping coefficient C = 0.25, 


and natural frequency w, = 10 radian/sec. 


a°a. dé. : 
+ 2cyus =z— +w 8 = 0 
ats n i; ar 





dt 


Changing to the required form (compared to Eq. (2-12) ) 








Wi) 
1 29% gees 7 
+. +42 —249 =0 
2 2 eh) Cie O 
a) atc n 
nN 
2 
1. Se So Sonne =" 0 
100 atc 20. Gt O 
Foe — = 0.1 
n 


Sa: 


e( 2) = - 2 ee 2a 
dt Ww dt 20. dit 20 


from which the equation for the Zero-slope isocline is 


obtained. 


(See Figs: (2-4) and (2-5)). 
3° "Discussion 

a. Any system which has a differential equation 
‘like Eq. (2-12) can use Liénard's method to construct a 
phase-plane direction field. After the direction field is 
constructed, one can draw any phase trajectory from a given 
initial’ condition on the phase-plane. 

b. As compared with the isocline method in Sec- 
tion A. If the differential equation can be manipulated 
to (2-12) form, a lot of labor may be saved, and one gets 
the same result. It seems better than the isocline method 
in this particular problem. 

c. If C is not a constant, then the zero-slope 
1socline will not be a straight line. 

d. Liénard's method (even for drawing an indivi- 
dual.trajectory) is not applicable to all second-order 


non-linear differential equations. 
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Fig. 2-4. Liénard's ‘construction 
of the phase-plane 


direction field. 
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Fig. 2-5. Direction field obtained 
by Liénard's method. 
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Cy DELTA METHOD 
1. Method Description 
Given 
¥ + G(x,x,t) = 0 (2-13) 
where (x, X, t) may be given analytically or graphically 
and may be a very complicated function. 
Rewrite Eq. (2-13) to the standa ra delta form. This is 


done by adding and subtracting the term p-x from Eq. (2-13). 
= ° 2 2 
x + (G(x, Xneet)eome Pox] actepsx = 0 


Let 


p65 a (G(x, X, Ey p°x] 


so that the linearity difference § becomes 


8 = +5 [e(x, & t) - pox] (2-14) 
p 


Eq. (2-14) will make (2-13) appear in the standard 4 form 


# + p*(x + 6) = 0 (2-15) 
The geometric significance of (2-15) in a phase-plane 


representation is seen by introducing the phase-plane co- 


ordinates 


oe and 


tS [xe 
ll 
< 


Then 


GWE pe dv .- Zeer 
Pat det ax 


* 
lI 
| 
ee 
II 
10: 
< 
| 
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The standard second-order equation (2-15) then becomes 


a first-order differential equation in x, v, and 4, namely 


Vi ev Gi Se 6 = () (2-16) 
avy _ he | a 
ax V (2-17) 


Fig. 2-6 shows in this case that the slope of the line. 
QP will be == ; consequently the equation of the normal 
to QP will be identical with Eq. (2-17). The instantaneous 
phase trajectory is therefore a circular arc passing through 
P and with its center at x = -6 or at Q. 


The next problem is therefore to relate the circular- 


arc segment beginning at P to the finite time interval At 


Since 

xX 

v= = 
p 

then 

CE cay 

at P 

Jon (2-18) 
pv 


Fig. 2-7. Let the average radius of curvature of the step 


be given by QP = 0. 


25 


x 


Normal to OP 





Flg. 2-O eee exact relations p. 





Fig. 2-7. Stepwise relationship 


‘fOr “Cons tale ap. 
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Then 


ds 
0 


eRe + dv 
ie hice ge 


2 
1 + (=) ax 


| (ea 
Page (ey 


—— 
— 


ox 
V 


dt = a qé (2-19) 
p 


Consequently, an integration will give 
tE- t - +(6-¢ ) 
p O 


pac = GAC (2-20) 


Eq. (2-20) shows, in regard to the angle subtended by the 
circular arc described by the phase point during the finite 
time interval At. 

A graphical integration of the differential equation 
in the standard delta form can then be effected by a 
successive construction of circular arc segments whose 
centers are located on the x-axis at various step value of 


x = -6& (see Fig. 2-8a). 


on 





(a) Change with positive 4. (b) Change with negative 4. 


Fig. 2-8a. 
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If x had been taken horizontal and positive to the 
right with v positive vertically, positive time would have 
been represented by a clockwise angular variation. (see 
Fig. 2-8b). 

The definition of 6 given by Eq. (2-14) transforms the 
original Eq. (2-13) into the standard delta form. Curves 
representing § in terms of the phase coordinates x and/or 
will be superimposed on the phase-plane so as to permit 
choice of a suitable average value of 6 for the construction 
of the circular-arc segments making up the phase trajectory. 

2. Example (I1I-B-2), Mass on Non-linear Spring 

The motion of a constant mass attached to a non- 


linear spring which becomes relatively stiffer with increas- 


ed deflection is described by an equation of the form 


Z 
Se 4 25(1 + 0.1 OS ae 
at 


Find a phase-plane solution curve for this equation by 


6 method, with the initial condition that at t = 0, x = 3 


Gr 
and ian Ox 


Comparing to the standard delta-form (Eq. (2-15) ) 


p- (x+6) = 25(1+0.1 ea ee 
as 
6 = Typuiises 


Sc 


15 





(a) Change with positive 6. 


Bag. 
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(b) Change with negative 6. 


2-8b. 


First draw the § curve. Phase trajectory construction 
proceeds as shown in the figure. The initial point is 
x(0O) = 3 and v(O) = 0, as the first step. It is assumed 


(1) _ 


that x "decreases to the vale 2.9 or, as -0.2 as shown. 
The average value of x during this interval is x__ (1) = 239% 


for which the average value of 6 can be read directly as 


we Gb 
V 


4 = 2.4 approximately. Thus the center of the frrst 


circular arc is located at point x = “6, B= -~24 and v=0; 
the radius is R = 5.4 and with this radius an arc is drawn 
to the point where x = 2.8, indicated by the small dot 
marked x(0) + Ax (1) in the figure. By continuing in this 
manner, the entire solution curve can be built up as a se- 
quence of circular arcs. (Since the original equation in- 
volves no damping, the solution is a closed curve repre- 
senting a periodic oscillation. (See Fig. 2-9). 
3. Discussion 
a. The delta method is limited to construction of 
phase trajectories for second-order differential equations. 
b. In starting a solution by the isocline method, 
the entire plane must be fitted with line segments fixing 
the slope of a solution curve. If only a single solution 
curve is needed, only a few of these line segments are 
actually put to use. Considerable simplification would 
result if only information related directly to the desired 
Eee von curve were used. 
The delta (6) method leads more directly to the 


desired solution. 
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Phase-plane diagram constructed by 4 method. 
The curve for non-linear function 4(x) and 
the final solution curve are shown. Dots 
indicate junction of successive arcs used 


iO Ieee uo. SOL UeLoOn. 
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c. The construction of phase trajectories will be 
carried out for various magnitudes of steps. So the accura- 
cy may depend on different amounts of constructional work. 

4. Notes 

a. 6 may be either positive or negative. The 
following 6 curves for five systems may help to understand 
how plots of curves may be superposed on the phase-plane. 


(1) The physical pendulum (See Fig. 2-10) 
a E. 
K + p sin’x = 


vy = + [6(x)-p*x] = sin x - x 
p 


(2) A system with restoration proportional to 


the cube of the displacement (See Fig. 2-11) 


x + 8B e = 0 


(3) A linear restoration system with quadratic 


velocity damping (See Fig. 2-12) 


x +q 1k |x 4+ p*x =0 





x=ql= 


[xe 


§ -+ ease -— p*x]=-% |x 
Pp Pp 
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Eage 


Bale 


2-10. 


2-11. 





The operative displacement 64 for the 
physical pendulum plotted in the phase- 


plane. Centers Op and O. correspond 


Q 
to points P and Q. 





The operative displacement 6 for a 
system with restoration proportional 
to the cube of the displacement. 


Centers Op and O. correspond to 


Q 
points P and Q. 


Ad. 


(4) A system with positive linear and with 


negative quadratic restoration in x (See Fig. 2-13) 
x + px -alxI|x=0 


§ = +5 (p*x St |) | p°x)= - = IS | x 
p p 


(5) A system acted upon by nth-power velocity 
damping and restored by first and mth-power displacement re- 


storation (See Fig. 2-14) 


% 4a |3271 | oe px + Al mt lx = 0 





= % |X) 2 1x , By m-l 
= Saale) i) an 

Pp 
ea 


(b) The 6 curve of a system indicates the system's 
type and degree of non-linearity. Thus, if a system is 
linear, 6 will be zero unless an artificial parameter has 


been introduced deliberately. 


1S SZEGO METHOD 


1. Method Description 


Consider the second-order system: 


7 
(2-21) 
£(x,y) 


hie 
ll 
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Fig. 2-12. The operative displacement 4 for 
a system with quadratic velocity 
damping. Centers O, and O 


E Q 
correspond to points P and Q. 
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Fig. 2-13. The operative displacement 6 for 
a system with positive linear and 
negative quadratoc restoration in 
displacement. Center O, corre- 


P 
sponds to point P. 
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Fig. 2-14. The two operative displacements 


6) and 6, for a system with 
damping proportional to the nth 
power of the velocity and with 
restoration proportional to the 
first and the mth power of the 
displacement. Center O, corre- 


Q 
sponds to point Q. 
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It 1s desired to identify the slope of the trajectory 
of this system at every point of the phase space. We con- 
sider the topographical system; that is a system of non- 
intersecting differential curves, which are symmetric with 


respect to the origin. 


v = $(x,y) (2-22) 


In connection with the function (2-22), construct the 


curves defined by the equation 


dv _ ov dx, ov ay _ 

dt ex at ay at ° (2-23) 
and Eq. (2-21). 
These curves we calilled GComtact curves. Substituting acm 


(2-21) into (2-23) 





OV day 

oo Cl enn R58, 

sv cl amet y (2 
Oy dt 


At every point of the curve (2-24), the corresponding tra- 
jectory of the system (2-21) has the same slope as the curve 
of the family (2-22) at the same point; furthermore, the 
curves, representing in the phase-plane the solutions with 


odd-order of multiplicity of Eq. (2-23) divide the phase- 


av 
at 


Peary, =e has the same sign on both sides of the curves, 


plane into regions in which has fixed sign. On the con- 


representing the solution with even order of multiplicity 
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of Eq. (2-23). Of course, in the latter case the slope of 
the trajectories at every point on these curves is the same 
as the slope of the family (2-22). 

In regions where x is negative, the representative 
point of the system (2-21) describes a trajectory in the 
phase-plane which crosses the curves of the family curves 
(2-22) in the direction of decreasing v. 

In the regions in which we 1s positive the trajec- 
tories cross the curves of the family (2-22) in the direc- 
tion of increasing v. 

Choosing different families of curves (2-22) and plot- 
=r = 0 in the phase-plane, 


provides enough information to sketch the trajectories. 


ting the corresponding curves 


If the scalar function, v = ¢(x,y) changes sign in the 
phase-plane, i.e., if the equation v = ¢(x,y) has some real 
non-zero solution with odd-order of multiplicity, then the 


curves representing these solutions must also be considered 


dv 
dt * 


This must be done in order to take into account the varia- 


as boundaries between regions with different sign of 


dv 
ae 


the curves of the topographical system (2-22). Notice that 


EPonein sign of due to variation in the sign of 


the trajectories of the system (2-21) on the curves v = 0 
have not the same slope as curves of the family (2-22). 


Consider the particular topographical systems: 


vo oF (2-25) 
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EOL Which 


av 
and 
Vy ay, (2-27) 


Ler Wire 


ee oe (2-28) 


The subsequent families of curves are now considered as 


topographical systems. 


V, =x +y (2-29) 
— x? = y" (2-30) 
MAT re (2-31) 
Using oe 
when Vee QO 
sh a nO) (2-32) 
and from 
ae A fs gE gt Sm A 
dt ox Gt ay Glits dt 
then 
fcr = 0) (2-33) 
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from 


then 


and 


Por 


Also 


But 


considering Vy > & + y 
xX = y = 0 and 
= 
dt 
dv 
ee ax dy 
a. - de “a ae 
= 2%y + 2y#&e,y) = 0 


oS 
7 = x2 se 
3. y 
Me im O pO) x = y 
and xX = -y 
Eo) 
dt : 

Nie Sa ae ay 
dt ax Gt od dt 
” Co dy 
= Set | ee 
= 2xy - 2yf(x,y) 

= 0 


ahh 


TE is seen that v 


2 


= 0 


(2-34) 


(2-35) 


Hence 


and 


Ns 


Panaliy , eonsideriang 


= £ (sey) 


il 
@, 


4 = *Y, gives 


FOr x = O and y = 0 


and 
A 
dt : 
But 
Wv4 _ 9%4 ax , 2% ay 
dt ox eke oy dat 
ks 
=y sy + xX £(x7y) 
= 0 
iene 
ye = -xf(x,y) 
2. Example (II-D-2) 


Consider Van der Pol's equation: 


<= 


II 


y 2 Lees - x 


De 


(2-36) 


(2-37) 


(2-38) 


(2-39) 


| ate = x 
(a) ae 0 for Var 5 
1 -x 
dv 
2 _ A 
(b) aes 0 feu y =0 
MoS ool 
dv 
There iS no variation in the sigm of at between the two 
sides of the curve y" = 0. 
dv ., 
(c) a = 0 fOr y = 0 
——— 
1 - 7 
dv 
(d) zc = O FOr y = ay NE: 21) ah(x2-1)24 4] -1) a 4] 


The next step is to plot the curves (a),(b),(c) and 
(d), (see Figs. 2-15, 2-16, 2-17 and 2-18). 

In order to use the method, it is necessary to super- 
impose the curves drawn in Figs. 2-15, 2-16, 2-17, and 2-18 
as is shown in Fig. 2-19. 

dv . 

The best procedure now is to check the sign of ae 
(v = 1,2,3,4) numerically at one point of the phase-plane 
and utilize the result for identifying the sign of ea 





(i = 1,2,3,4) in the different regions in which curves 
hv. 
aaa =O and v = 0(i = 1,2,3,4) divide the phase-plane. 
Since 

Oi V5 < + © 

=ole < 

c Va + © 

dv5 dv. 
It Zs necessary, before checking the sign of ae? aimee 
V 

and aa at one point of the phase-plane, to identify the 


Veatever dt that pont. 


Syileigy 1onm Vir V3 A 
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generate regions 
M,N,f,O,R and@a 
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dv 
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aD 


Ze 
ae. S 0 ~SFOr 
y* = 0 
re Sl 


generate regions 


o,B and v 


V3 = Oo5 for tx 
dv 5 
ae O for y 


generate regions 
AVE Ce, :, F,G, oH, 


I,J, and K 








aA = 0 
[Orme = Of y= ao 
g h 
ae, 
dt 
+ 
f oe i) = | (x*-1) w/6e-1 +4 
A +[ 2 
= 0) x 
e : dt generate regions 
a,bac,d,e,£,9, anc 
d C 
AVG. V4 =O 
dt? 
dv 
Pcs 2k. IG es Cee O and or a OF 





Fig. 2-19. Example (II-D-2) 





Some information about the slope of the trajectories 
at every point of the phase-plane is thus provided. For 
example, consider the point P(-3; 1.5) in Fig. 2-20), 


Giy 








from = <0 , requires that the trajectory move into 
y the region A'F'D'B' 
dv . 
from Ae >0O , the trajectory moves from P into the 
- region E'C'A'F' 
dv, 
from ac <0 , requares that the trajegitory move anto 





the region C'A'F'D', 

These three conditions, taken simultaneously, define 
the angle A'PF' in which the trajectory which originates 
at P must lie. 


Since the slope of the trajectories is known exactly 
dv. 
at the points on the curves =a = 0, it is possible to 


draw the trajectories starting from P with good accuracy. 
(See Fig. 2-19) 


Starting for example from point B(-3,1.5) contained 
av. 


in the region of fixed sign of ae » namely T, @, I and G. 


(See Fig. 2-19) 
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vz 
=O o4=0 


Te 


JV, 
570 


“ © O | 
bY O > 








av 
5 
dt alt 
Pp 
dv 
4 
= < 
at o 
P 
NOTE: 
dv, 
—_ < o 
AL O ism T, OO, M and S 
dv 
—_—— YE 
dt O in N and R 
dv. 
———  <— , 
dt O iain @and ¥ 
dv, 
ae O in 
dv, 
—— <0 in BS Dr, H, J atid 
ert . 
dv. 
ay > O in @& C, EE, G, L ame K 
ie 
dv 
a7? in tt, d, f£-anavk 
dv, 
—S  < ; 
at O ain a,c, @ and g 


3. Discussion 
a. This method is a simple and useful method for 
drawing trajectories. Also it is essentially a generaliza- 


tion of the isocline method. 


a 


bs Accuracy Gan be improved by considering further 
topographical systems. 

c. It is worthwhile to notice that the isocline 
method is a particular case of this method. And it corre~ 


sponds to the particular topographical system: 


Vv = a,x + a5y 


NOTE: 


N,N Ea Ly 
dt ~ 35 dt sy - dt 





= ox dy 
Sac fo 3 
= 0 
dy ” 
(Agee So eae 
dx Cc Dal a5 
dt 


Bis MURTHY'S NEW APPROACH TO THE PLOTTING OF PHASE-PLANE 
TRAJECTORIES 
1. Method-Description 
This method is a new and simple approach for the 
isocline plotting of second-order time-implicit non-linear 
differential equations given by x = F(x,x'). 


a. Isoclines of the equation x"+f(x')+gq(x) = 0: 


x" + £(x') + g(x) = 0 (2-40a) 
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Substituting Mx’ = x" anmboe(2-244) 


dx ! 
where M = dx * gives 


Mx" + fx") + g{x) = 0 (2-40b) 


(1) Plot from the given differential equation 
(2-40a) the curves of x" vs x' on the x"x' plane for con- 
stant (different) values of x. This is a simple matter 
Since it amounts to shifting the curves of x" = £(x') 
vertically for different values of x. 

(2) Draw a number of straight lines x" = Mx' 
through the origin (corresponding to the different values 
pf M) on the x"x' plane; this can be done quickly. 

(3) With any particular straight line in 
view, say for M = M,. determine all possible points of 
intersections of this straight line with the curves obtain- 
ed from (a). Thus a number of values of x and x' can be 


obtained corresponding to M=M which when plotted on the 


af? 
x'x plane and joined together represent the isocline for 


5° M, ~ oe Witt 


a fair number of isocline curves are obtained on the phase 


M=M,- The procedure is repeated for M=M 


plane. The respective slope lines are marked on these iso- 
clines. Joining these slope lines gives the different tra- 
jectories or the entire phase portrait for Eq. (2-40a). 


(b) Isoclines for the equation 


x" + £(x')h(x) + g(x) = 0 (2-41) 


ol 


The above procedure can still be applied except that the 
curves of x" = £(x') are now multiplied by a scale factor 
[due to the h(x) term] corresponding to different values 
Oi. x. 
(c) Notes 
(1) This method is more convenient in deter- 
mining the solution curves for the following differential 


equations: 


x"+ x'h(x) + g(x) = 0 
5c es ao) — a0 (2-42) 
x" + xx! cna — 0 


when compared with existing methods. 
(2) Thisemethod can also,be,yused for the 


following nonlinear differential equations: 


O 


h(x)x" + R(x) £(x') + g(x) 


(2-43) 


ll 
S 


tx )xe re KUX) £(x") teg (x) 
Salus: Filscemesc.')) 


(3) The Lienard method is a step-by-step con- 
struction for drawing any particular trajectory correspond- 
ing to a given set of initial conditions. The present 
method determines the entire phase portrait indicating the 
Dold Victor ecil, initia a leeondlegons: 

(4) Variations of the nonlinear restoring- 


force term g(x) can easily be studied. 
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(5) The method is independent of the nature 
of the functions f,g,h,k, i.e., whether logarithmic, alge- 
braic, or transcendental. 

(6) For nonlinear differential equations of 
the type x" = f(x), the x"-x' curves become horizontal 
straight lines for constant values of x. The intersections 
of these straight lines with x" = Mx' lines give the points 
on the phase plane for any chosen isocline slope. 

(7) The isoclines corresponding to M = O and 
M= © can generally be plotted quickly. 

(8) Sometimes an x"-x' curve (for a constant 
x) and an x" = Mx' line may run so close to one another 
that the isocline curve on the phase plane is nearly ver- 


tical in that range of x'. 


2. Example 
Consider a system with nonlinear damping and re- 
storing forces described by x" + ie “5 x? = 0. 


(a) Draw x"-x' curves (see Fig. 2-22). 

(b) Draw Mx' = x' corresponding to various 
values of M (a number of straight lines through the origin. 
(See Fig. 2-22) 


(c) Draw isoclines and phase trajectories for 
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Considering the line for, say M = + 5, the points of in- 


tersections are determined as 


c= 0 jee ee) 
(x= 2 , x= =) 
(x = /3°, x' = -7) 
a= 9 = 
@=/6 , x" = -2) 
(x = /6 , x' = -2.5) 


These points are marked on the phase-plane, joined and the 
Slope lines are marked on it (Fig. 2-23). Similarly, dif- 
ferent isoclines are obtained for different values of M. 
The trajectories are then drawn as shown in Fig. 2-23. 

3. Discussion 

(a) This is a simple, all graphical, useful and 

general method to obtain the phase-plane trajectories cor- 
responding to all pertinent initial conditions of second- 
order autonomous systems (which are not amenable to the 
usual isocline method) by making use of the isocline 
method. It reduces the labor of calculating isoclines for 
many non-linear differential equations and thus makes the 


i1socline method a much more useful todl. 
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be 





-3 


Fig. 2-23. Isoclines and phase trajec- 


tories for x" + xn? “+ — Oe 
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ine THE ACCELERATION-PLANE METHOD 


1. Method-Description 


It is based on the fundamental equation 


dav 
ax 


’ (2-44) 


<1 


where w = w(t, v, x) or the acceleration is a function of 
t, v and x. For the sake of simplicity, assume that a con- 


stant driving force is applied. 


Thus dv _ w(v,x 


a = (2-45) 


The method calls for preliminary plots of w vs x for 
constant values of v. This auxiliary w-x plot for constant 
values of v iS not a trajectory. It characterizes the sys- 
tem just as the original differential equation character- 
izes the physical condition of the system. 

This method is applicable to both linear and non-linear 
problems. In the w-x plot, though constant driving func- 
tions can be included, the time-varying forcing or driving 


functions must be added for different values of time. 


2. Example x + OF 5° eee — oe 
Forw= z' = 0, then x” === ,»and tree” —- =x 


plot or w-x plot is a straight line with a slope -l . 
For v = x' = 1, then x" = -x -0.5, so this straight line is 
displaced from the first line by -0.5. 

In Fig. 2-24(a) six points of construction are shown 
for the v-x trajectory, starting from PS where t = 0, 


<i 2.5 ane v = 90. 
Oo Oo 
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Taking xX, as the center, and starting from Wi, an arc 
can be drawn in the clockwise direction so that the arc 


strikes Ry on the horizontal axis. From Ry draw a line to 


Pi where x = x and v = Vii the perpendicular line to 


PiR) gives the correct slope for the v-x trajectory at P,- 


If the perpendicular line gives a slope that agrees fairly 


well with the assumed slope, point P, can be considered as 


Hl 
the desired point which falls on the v-x trajectory. 
WwW 
Slope of perpendicular line = —_ ‘ 
I 


If the assumed slope does not agree with the slope of 
the perpendicular line obtained from the above construction, 
a new slope should be assumed and the above procedure re- 
peated until the assumed slope agrees with the slope of the 


perpendicular line (see Fig. 2-25). On Fig. 2-24(a), P 


1 
has the walue’x = xy ana vy = Vy and Vi is chosen as -0.5. 
Thus if one prolongs the vertical line from xy downward 
through Pj, this line should intersect with the v = -(%) 


line in w-x plot (shown in dot lines) at a point W,- By 
connecting the points Wor Wye Wo eee etc., the w-x trajec- 
tory is constructed at the same time (see Fig. 2-24(b)). 
3. Discussion 
This method is applicable to both linear and non- 
linear problems. One draws the w-x trajectory first, then 
gets the v-x phase-plane trajectory. Each time only one 


trajectory is obtained, so this is a useful method. 
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W-X 
EQUI-VELOCITY 
LINES 


Fig. 2-25. Acceleration-plane method. 
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oe DEEKSHATULU'S METHOD 
This method suggests the use of simple transformations 
for second-order non-linear differential equations to effect 
rapid plotting of the phase-plane trajectories. 
1. Method-Description 
a. Transformations for second-order autonomous 
systems 
A general second-order non-linear autonomous 


system is governed by a differential equation of the form 


eee (oe), (2-46) 
Let 
e dag 
Y= at 
then 
2 


dt 
_ ax 
ED meats 
Change Eq. (2-46) to 
Pie sa fel (cases) (2-47) 


To draw the isocline for an assumed value of M, equal 
to M,, say, x has to be found as a function of x and M, and 
this in many cases is not possible; thus the only attempt 
made is to find the points of intersections of such an 


1socline with a given set of straight lines or with simple 


fA 


standard curves. Consequently, on the phase-plane, if one 
considers the intersections of the isocline for slope M=M 
and the set of straight lines x = Kx(K = +1, K = +2, 


Eq. (2-47) becomes 


M, Kx =i Kx) 


ONG 


M, = F(x,K) (2-48) 


For a specified K, the solution of Eq. (2-48) (algebraic 
equation with constant coefficients) for real roots of x 
gives the points of intersection. But it is easier to find 
the values of the slope M (for various values of x) on the 
lines x = Kx than it is to try to obtaagh the isocline eum. 
for a given value of the slope. Even under this transforma- 
tion x = Kx, every quadrant of the phase-plane is uniquely 
defined. It is advisable to find any symmetry present in 
the phase portrait for a given differential equation since, 
if it exists, it reduces labor and is tested as follows: 


(1) Symmetry about the x-axis 


Tf for K = positive (negative) Ky and 


Ds positive (negative) x,, the slope M = My. then for 


K 


negative (positive) K and x = positive (negative) Xi 


i ba 

the slope M = -M 
(2) Symmetry about the x-axis 

If for K = positive (negative) K, and x = 


positive (negative)x,, the slope M=M then for K = nega- 


a 
tive (positive)K,, and x = negative (positive)x,, the slope 
M = -M>5. 
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The procedure for determining M, on the phase-plane 
for any chosen value Ky and x, Can be viewed in the phase 
space (with x, x, and X as coordinate axes). Determine the 
curve of intersection of the surface defined by Eq. (2-46) 
Vwithethes planes = Kx, then determine the simple point of 
intersection P of this curve with the plane x = X)° The 
slope of the plane passing through the x-axis and the point 


P with x-x plane determines the value of M Similarly, 


1° 
other values of M are obtained by taking different values 
of x and the procedure is repeated for each value of K. 

Other transformations such as x = ae x = Kx>, x or 
x = K, etc., can also be considered. 

This procedure is general, systematic, and comparative- 
ly simple when compared with the method of determining slope 
lines arbitrarily at points for the phase portrait construc- 
fren). 

b. New Planes for Certain Second-order Time- 

Varying Systems. 

Time-varying elements occurring in physical 
systems; for example, the variable mass of a guided missile 
can be though of as function of time since its weight de- 
creases as the fuel is consumed. 

On the phase-plane, solution of this kind of problem 
makes use of step-by-step procedures by considering simul- 


taneously the x-xX and x-t (solution) planes. “~Proressex 


Deekshatulu suggests some new planes whose coordinate axes 


te 


are explicit functions of time and which effect an easy 
and rapid solution of at least certain classes of second- 
order non-linear time-varying systems not easily amenable 
otherwise by the phase-plane methods. 
c. Suggestions for Solving Certain Higher-order 
systems 


Assuming non-linear transformations of the 


EOaan 


@ 


f=et+af(e) or n= 6 + bf(e)+ aé (2-49) 


third and fourth-order time-varying systems described by 


the following equations: 


c-E + £ + £(6) 


= OQ 

ee ee E) =O 
(2-30) 

Sa eee ne (Cah) 0. 

2th +n + E(n) = 0 


can first be solved by the method described and then be 
solved as a first-order or second-order nonautonomous sys- 
tem. Note that the equation of the form é + af(e) = h(t) 
can be solved on the e-t (solution) plane. 
2. Examples 
(Example 1) Consider a non-linear system with dis- 


Sipation proportional to xf, given by x + hx? +x = 0 
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Let 


dx 


dt 
then 
Z 
dm. My _ dy ax .., Wax 
2 dt Gx) “Git Ww Gt 
dt 
Thus 
x + te + x 
™ a2 
=Srix + Ax + x = 0 
ay a 
og hx +Mx + x = 0 


Both roots of this equation are valid. On the other hand, 


if one substitutes x = Kx this equation becomes 


Ma fe oe 6 
or 


eae 


M = K 


which is simpler to plot for any given values of h and K. 
Since the phase trajectories are found to be symmetrical 
about the x-axis, only the upper half of the portrait is 


shown in Fig. 2-26 for the case h = l. 
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A 
il 
pa 
Me 
lI 
*% 


at 


M = 2x - 0.5 


etc. 


(Example 2) Let it be required to find the behavior 


of the following time-varying system: 


aoe + tx + = = 0 


for all pertinent initial conditions (with t # 0). 
The task would be enormous if it were to be attacked 
by the phase-plane method. Using the tx-x plane the dif- 


ferential equation reduces to 


x? 
oo = = sue 
Note: 
Let 
f= y 
then 
om ._ Gy Ce _ za 
dt a oe Mx xX + tx 
: ‘ 3 
Mtx = tx + t X = =x 
or 
x? 
x = y= = “va 


v7 


Taking different values for M the isoclines (enlarged 


or reduced scale plot of the cubic curve) -are drawn and the 


respective slopes marked on them with their. direction (see 


method, 


°T> 


EV ciane 2 = 27 ee - 
3.) Discussion A 
a 


a. Simple substitutions like x = Kx, Kx", 


~ LD, 2 Z « . ° e " « ‘é “ = N : 
x +x = y ain the isocline equations of ‘second-order sys- 


tems Simplify plotting of the phase trajectories to a great 


rN 
J: 


extent. The method is a useful supplement of the isocline 


oie si eehs aH ee ; Es ae oo ees et =. i Dt dm “4 
DONtI oF Beir row: Sioa: Sha {Sc etaoneue 


b. Sinkgathers PEs POS Te7y metbege to, evaluate 

quickly the bahavior of second-order non-linear time-vary- 
¢ . 

ing systems for various pertinént initial conditions and/or 

for’ ‘dif feréntynon-linear characteristics;,;Deekshatulu 

suggests), an=general time-varying coordinate,planes, the 


° SZ -. ‘ 
use Of; tx-x--ands tk sx planes;:fox; obtaining, easily the solu- 


tion of certain classes of nonlinear. time-varying, systems. 


H. PELL'S METHOD = ey 


1. Method-Description 


" 


Quite generally, second-order nonlinear autonomous 


equations of the form - 
2 
clic 


may be written as 


a“x ae Le. _ 
f= + 6 (=) + E(x) = 0, (2-3ay) 
2 che 

dt 10 

c 

~y | 

a ate j we} 
M | 
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Let 


_ Gx 
Y ~ GE 
then 
2 
da x dy Tae dx 
dt Get ax Gite 
Thus 
d°x ass 
—>5 + 6 Ge ate 
dt 
. gy, dx Ox 4 
ae edie b(5,) + £(x) O 
or 
ax 
dx ox 
dt 
a cada Ee (2-32) 
VY = 


Now, as shown in Fig. 2-28, plot -é(y) versus y and 
-f(x) versus x. Given the initial conditions x(0) and y(0), 
the construction of a segment of the trajectory through the 
initial condition proceeds as follows: 

a. Move vertically from IC and find f[x(o)] ; 
lay this off on the x-axis from x(0) by use of a compass or 
a 45° triangle. This establishes point A. 

b. Draw line AB, where OB = y(0). 

c. Find -¢[x(0)] by moving horizontally and draw 


CD parallel to AB. 
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d. Now the hypotenuse (D) (IC) of a right triangle 
whose base is the distance DE = ¢[y(0)] + £[{x(0)] and whose 
altitude is the distance (E) (IC) = y(0) is established. 


The slope of the hypotenuse is 


Bl cy _ y (0 
DE ~~ ‘ly (oy) + £tx(0) (2-3 


e. The negative reciprocal of Eq. (2-33) is the 
slope required in Eq. (2-22); . thus a line perpendicular to 
the hypotenuse, passing through the initial condition forms 
a segment of the required trajectory. The usual precautions 
must be observed as to We 1evrettn of the extrapolated seg- 
ment. 

Note: The major advantages of this method are: 

(1) Its simplicity. Only two triangles are 
needed to carry out the solution once the two functions are 
plotted. 

(2) It eliminates the trial and error of the 
delta method. 


2. Example 


Consider the equation (the hydraulic speed control 


system) 
ace de 2 Zee Ke 
—=—-+e,—-+ ua “(ltav°e’ )e + == 0 (A) 
2 liga O a: 
dt 2 
where 
c, = viscous damping coefficient 


om 


c, = inertial coefficient related to volume of nozzle 


Z 
chambers 
K = gain of amplifier 
w= natuisal frequency of vibration of flapper tem 


small amplitudes. 
The equivalent flapper spring is assumed to be linear. 
Eq. (A) may be rewritten as 


d“e de 


2 2 ae 
late? Ww, (1+b e-)e = 0 (B) 


Let T= Wt to change time scale; thus Eq.(B) becomes 


wipe yeas GS 
2+ @ at 


; 


Glen eave =k) (C) 
Git 


kK 


Choose the following more or less arbitrary set of 


coefficients designed to emphasize the construction tech- 


nique: 
c, = 164 W = 377 
Il O 
2 
Cy = O50: a = 0.858 
K = 1000 


on a consistent set of units. Then 


2 a Sia 
w* = w” + a 242000 
2 
b’ = 2. = 0.58 
7, 
a’ = On 
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Let 


2 SSS 
_ de 
Y ~ at 


Eq. (C) becomes 


dy _ -0.33y -(1+ 0.58x“)x 
dx V 
éd(y) and f(x) may be identified and plotted and the 
construction completed as shown in Fig. 2-29. 
3. Discussion 
This method is satisfactory: 
a. If only a single solution curve is needed. 
b. Only a few of these line segments are actually 


De CO use, 


ine SOMAYAJULU'S METHOD 

In this method a graphical procedure is suggested for 
drawing phase-plane trajectories making use of tangent 
slopes instead of normals (i.e., the negative reciprocal of 
tangent slope values) as is done now in the Liénard con- 
Steric bon . 

1. Method-Description 

A straightforward construction is presented for 

more general cases to which the isocline method is not 


applicable and construction of individual trajectories is 


required. 


2),5: 






Cts 0), 4(0)) 
| 
| 
iF 

fx (0) 












vi * 





f(x) 
Fig. 2-28. Pell's Method of phase.trajectory 


COonusSeructton ; 


CB=AD 





Fig. 2-29. Example construction using Pell's 


method for the hard-spring system. 


84 


a. Let the slope of a second-order system be 





given by 
_ GX _x 
aes 


£, (x) 4 £.. (x) 
g, (x) + 95 (%) Foe 
Assume ti 


and of the odd characteristic type; these curves are first 


fo, Gy: and G5 are all continuous, single-valued 


drawn on the phase-plane. 

Referring to Fig. 2-30, the initial point P(x, +X.) is 
marked. Make PQ = £,(x)), QR = £,(x,)- Similarly the 
length RST is made equal to J, (x,) + 35 (X))- The angle PRT 
1S a right angle. Note that the horizontal line RST is 
drawn to the left of PR since the slope is Known to be po- 
Sitive from the nature of the functions {1 fo. Sy: and Sg. 
and the initial conditions oe and Kos A small length of the 
straight line joining P and T gives the required slope at P. 
The field direction is marked on the basis that with x posi- 
tive x should increase. It is seen in the procedure that 
it is necessary to first determine the direction of the 
normal to the phase trajectory. 


b. Consider equation 
x + £(x) h(x) + g(x) = 0 (2-35) 


This equation is considered to be of practical impor- 
tance because of the presence of nonlinear damping and 


nonlinear restoring forces. 
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Let 


dx 
Gt 


and 


rewrite the equation (2-35) in form (2-36) 


ee a (2-36) 
x 


The plots of f(x), h(x) and g(x) are shown in Fig. 2-31. 


Also a parabola x7 One x? (only one half will suffice) is 


drawn in addition to the above curves. 


° 


Let the initial point be at P(x, 


then: 
OS = XG » OoOdcH= Xo 
OR = h(x) ; obo = E(x)) 
and SA = 


g(x.) 

Make SM = QR. 

Draw a circle with BM as a diameter and let it intersect 
the x-axis at the points C, C. 


Then 


sc = JE (x) h(x)) 


With the help of the parabola the squared length of SC 


is obtained and is marked as PN. 


Make NL = g(x) 


and LK = x 
O 


So? 


“s—TTeqjep uoTqJONAASUOD “TE-z *6ta 





CX) 5 q 


<— +X 


(x) U 


(X) 6 


Og Voy 
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The segment line through P joining P and K represents the 
required tangent slope of the phase trajectory at P. The 
direction of motion is decided to be downwards. The new 
values xy and Xy at the end of this segment are determined 
and the procedure repeated. 

Note: (a) In both the above mentioned methods it is 
necessary to take equal scales on the x and x axes. 

(b) The accuracy in plotting is enhanced by 


considering small segments. 


2. Example 


with initial condition Xo = -4 and ie = 0. 


Rewrite x + 7 +x = 0 to the form 


- - x Bo oe 


M = 2 
x -X 
and draw 
= ae = 
fi = f., — i 
Crile 


on the phase-plane (see Fig. 2-32) (using Somayajulu's 
method as in I-l-a). 
3. Discussion 
This method is based on geometrical relationships, 


SO a precise step-by-step procedure must be followed. 


eo 





oT durex qq 


“Ge c 


*bTa 
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ra THE E-FUNCTION METHOD 
1. Method-Description 


q7 ax 
Se + £(= , x) + g(x) = h(t) 
ate dt 


(2-45) 


Let the Kinetic energy and potential energy be defined as 


eZ 


(4)x" and Q(x) 


where 
x 
Q(x) = | SMe.) aoe 
O 


Define the total energy as 


2 


E = 5(k)" + Q(x). 


Differentiating E with respect to time 


2: 


DE _ dx, d x , dQ(x) 
Siege Gla dt 


R=x - xX + Q(x) 


Differentiating E with respect to x 


Il 
re 
=e 
iQ 
x 


(2-46) 


(2-47) 


(2-48) 


(2-49) 


From Eq. (2-49) 


d x dE 
ap ee) 


Substituting this into Eq. (2-45) yields 


SE - g(x) + £(%, x) + g(x) = h(t), 
OE 
= = Ge) = 6 (Gc 2) (2-50) 


From Eq. (2-47), an E-x plot can be made for different 
values of x. These curves may be called equi-velocity 
curves, as xX = v, and they represent different kinetic 
energy levels. For x = 0, the E-x curve is simply E = Q(x) 


where Q(x) represents the potential energy, At x = x,, the 


ne 
difference between any two E-x curves, say, between the 

curves for v = vi and v = Vo 1s equal to the difference 
2 2 

Lao Ue 


The value of v = x is, from Eq. (2-47) 


x =v = + /2[E-Q(x) J. (2-51) 


(ye 


Substituting Bq. (2-51) Tnto (2-50), 


dE 


ae WCE) -£| "x, +/2TB-0(x) J | fen vs Oo 


(2-52) 


= h(t) -£| x, -/2tB-9(x) J | for v <0 


o2 


The E-function method is thus a method for construct- 
ing the E-x trajectory according to Eq. (2-50). See Fig. 


2-33, the slope of AB is equal to h(t)-£(x,v_.), where 


_ Waa: 
Vav 2 
2G Xo 
and x should be corrected to “az OeVen though Peas 


initially taken as equal to X4° 
2. Example: Van der Pol Equation 
g-(1 - x*)k +x = 0 
E-x plot (see Fig. 2-34). 
From E-x plot get the v(v = x) 
Note: As compared with the isoclines and trajectories 


of the Van der Pol equation obtained from the digital com- 


puter, errors exist in this method. (See Fig. 2-35) 


i THE SLOPE-LINE METHOD (See Fig. 2-36) 


Given an equation 


coh 
asc £ (x) 
or dy = £(x)dx, the incremental relations may be written 
Ay = fy 4% (2-53) 
ft, + f. 
where fiy = represents the average value of f(x) 


during the interval. 
Thus 
By See ee  E 
(2-54) 


lI 
S 
= 
~ 
sf 
NO 


a5 


Vevy 
¢ 
A 
/ 
/ 
8 | VeV, 
/ 
ARE ~~ V=0 
: K.E. Je 
| | Q(x) 
| PE 
Xy Xe O aa 
V+ 
Adsl cOenin-f (x, oe 


Fig. 2-33. The E-function method. 





Fig. 2-34. Example (II-J-2). 
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Let an angle 9 be chosen such that tan 6 = oe for 


a chosen value of Ax, then 


> 


-e 
T 


tan 8 = 


Ay. 
am (2-55) 
Starting from an ordinate Ei: make an angle 9 and draw 
a dotted line bc to intersect the horizontal line atc. 
Then ac = Ay. From c, draw another dotted line, making 
the same angle Gmwmen the vertical line erected from c, 
until another ordinate f, is intersected. As shown de = f 


2 2 
can be projected from the f(x) curve on the left of the 


figure. 


Then ce = Ay. and ae = Ay, corresponding to the in- 
Grement Ax shown on the left of: the figure. 

This method can be applied to the solution of simulta- 
neous first-order differential equations, linear or non- 


linear. For a first-order nonlinear equation 


GY + g(x) = h(x) (2-56) 


the incremental relations give 


Ay = [h(x) - g(y)_ JAx (2-57) 


2 


0) 
0) 
G 
A 
“ 
U 
O 
Uv) 
4 


and 


trajectories 





. (Example II-J-2) 


2-35 


1a eo 


a. | 





at 


A 


The slope-line method. 


2-36. 


Fig. 
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For example, in a problem of surge tank transients, the 


equations are of the following form: 


- Bay + £(u) (A) 
- Y= v7, y) - u (B) 


When the equations are in normalized form both the 
horizontal and vertical systems of slope lines can be 
drawn with the same slope, Eq.(A) represents acceleration. 


Eq.(B) is obtained from the condition of continuity. 


L. COMPARISON 

It is readily seen that some methods are easier than 
others but may be of limited usefulness, applying only to 
some particular problems. Some of them are useful but com- 
plicated, so that it is easy to make a mistake. It seems 
the most useful and general method is the isocline method, 
when one solves problems by hand labor. Deekshatulu's and 
Murthy's methods are useful supplements. The difficulty 
of the isocline method lies in the labor required to get 
the isoclines. So if we can use the digital computer to 
generate a program, call upon it to draw the isoclines and 
phase trajectories on the same phase-plane, then the iso- 
cline method is more useful and general than other methods. 

For time-varying systems, the Delta method and Pell's 
method are useful. Deeskshatulu's transformation method 


(by using a new plane) is a useful supplement. 
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III. THE APPLICATION OF PHASE-PLANE METHODS 


Phase-plane analysis of nonlinear system problems not 
only gives a convenient display for interpreting results, 
but also permits use of different techniques for shaping 
the phase trajectory to make it satisfy design requirements. 
So, it can be said that the phase plane has become a logic- 
al tool in the solution of problems. 

Generally, phase-plane analysis is limited to second- 
order systems. Since higher-order derivatives cannot be 
plotted on the phase plane, higher-order systems would not 
be completely defined. Of course, the method might be ex- 
tended from phase plane to phase spaces of n dimensions to 
handle nth-order systems, but then it loses some of its 
simplicity and convenience. 

In this chapter, the manipulations of the phase plane 
are discussed, then the uses of different graphical methods 
and techniques for phase-plane analysis of nonlinear system 
problems. All the problems are of second order or can be 


approximated by second-order differential equations. 


ie MANIPULATIONS ON THE PHASE PLANE 


1. Elementary Algebraic Manipulations 
a. Addition 
Addition of phase-plane trajectories is accom- 


plished ina fairly straightforward manner. Suppose we have 
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two quantities xy and Xoo and it is desired to find the 


phase-plane trajectory of xy + Xo. 


Let x and x represent phase-plane coordinates: 


therefore 
x, = xy + Xo 
; : ; (4-1) 
xy = xy + X5 


On the graph, simply add the x's to obtain the new x at a 
particular time and the x's to obtain the new x. The time 
for each plot must be known at each point so that addition 
takes place at points of isochronism. 

bere Mul ei pi reation: 


If the trajectories are x, and x, then the 


1 2 
products are: 
KX, =X) ° Xp (4-3) 
KX, = X)X> + X5X) (4-4) 


Again the addition and multiplication must take place 


at points of isochronism. 
c. Raising to a power: 


For nth power 
x = x (4-5) 


x, = nx x (4-6) 


he, 


The plot is taken point by point and the x coordinate 


1s raised to the nth power to obtain the new Xi the x co- 


ie power and multiplied by nx 


Ordinate Vee rarsed@rerethe x 
in ordes to obtain the new Xe 
d. Operation on the trajectory with a function: 


EE the trajectory for x is Known, then “ene 


trajectory for f(x) can be founded as follows: 


Pe £ (x) (4-7) 
idee dt 
(4-8) 


ax 


The new x is the value of f(x) and the new x is the 
slope of the f(x)-vs-x curve times x. 
2. Differentiation 
The trajectory x 1s Known and the trajectory of 
the derivative of x is desired. This has a practical appli- 
cation. Suppose a device or system output is to be display- 


ed as a phase-plane trajectory. This can be done by differ- 


ox 
dt 


is plotted against x. Unfortunately, it is very difficult 


entiating the output x, thus obtaining or x. Then x 


to differentiate properly; and alternate is to take the in- 


tegral of the trajectory, then differentiate the trajectory 


graphically to obtain the true x-x plot. 
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For illustration purposes, let 


u= [xde; (4-9) 
where 
x = du 
atin 
Then 
3 = KK 
dt 
But 
gx _ dx du _ dx 


dt dau dt dau * 


So in order to find the new x for each point x and u: 
a. Find the slope a : 
b. Multiply this. value by the value of x at that 
point . 
c. Plot the new x on the -u axis. (See Fig. 3-1) 
3. Transformation 
It is possible to transform the phase-plane tra- 
jectory by applying some transformation to x and x quanti- 


ties for some purpose. 


For example: 
ia £(X,. x4) 


x5 = G(X). x1) 


An example of this might be 


X, = X. sin x 
2 al i 


* 
Cd ry 





Fig. 3-l. Differentiation procedure. 


G2 





4. Methods of Shaping Phase-Plane Response 


There are several ways in which a phase-plane re- 
Sponse can be shaped to conform to design requirements. 

These are: 

(1) Choosing specific nonlinearities. 

(2) )GLimi tinge 

(3) Piecewise shaping. 
a. Choosing Specific Nonlinearities: 

For the following second-order differential 


equation: 


x + £(x, x)x + g(x, x) = 0 


various combinations f(x,x) and g(x,x), may be chosen to 
obtain a required response. 
bo) Limverng: 

The values of x and x may be restricted to 
certain values as shown in Fig. 3-2. The limiting may be 
extended to x also or even to some function of x and x. 

c. Piecewise Shaping: 


See Fig. 3-3. 


Br. COMMON PHYSICAL NONLINEARITIES 
1. Saturation: (the most common of all nonlinear 
phenomena) 


Consider the case of a saturating electronic ampli- 
fier driving a motor in an instrument servo. (See Fig. 


3-4(a) and (b)) 


GS 


Limur 





Fig. 3-2. Limiting of x on the phase-plane. 
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Fig. 3-3. Piecewise shaping of trajectory. 
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(a) Block diagram of a servo 


with saturation, 


Output. 


Linear 
approximation — 
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(b) Input vs output characteristics. 


Fig. 3-4. Common physical nonlinearities - Saturation, 


EOS 


The equation in the linear region is 


«KK K 
E+——"E5=0 


eo 1 
_ — 
= T NT 
m 


il 1 
For saturation operation: 
ee ae ee eG Bee 5 
7% 1 
m 
eee ee E2+E 
us 1 
m 
where 
15 
af Ee 
NT 1° 
m 


Using the isocline method, draw isoclines on the phase- 
plane E vs E. In the linear region the isoclines are radial 
straight lines (see Fig.3-4c), emanating from a focus at 

the origin. 

While in the saturation region the isoclines are hori- 
zontal, the transition from linear to nonlinear operation is 
indicated by a vertical line at E = E,: Dividing lines on 
the phase plane may be obtained from many nonlinear phenom- 


end. 
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Fig. 3-5. Dead zone - (a) Input vs output character- 


istics (b) Block diagram of a servo with 
dead zone (c) Isoclines, dividing lines 
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2. Dead Zone 
This arises in mechanisms which are spring-loaded 
to i: eentiee backlash and in many other devices which are 
insenSitive to small signals. Consider a servo with a dead 


zone adn the error transducer. 


E+ f+ =0 
Tm 
= Ey < E< + rT 
ee eS 
7 Nr 1 
le |> 5, 


Tn the dead zone, there is no control forcing-function, 
so in this case isoclines are straight lines with slope 
-pop 
Tm 
M= 2 
E; 
(See Fig. 3-5). 
3. Backlash 
Most mechanical linkages exhibit a certain amount 
of backlash. (See Fig. 3-6) 
Consider a simple positioning servo with backlash. 
When the backlash is taken up, the operation is linear. The 


eguation of the system is 


E + Bic iso + i o =.0 
n n 


1S. 
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Fig. 5=o7—™ Common physical nonlinearities = 
Backlash. 
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When zero velocity is reached and the motor reverses, the 
output shaft remains stationary until the backlash is taken 
up. Thus the error remains fixed at some value E,- 


The eguation of this motion is 


a 2tu 6 = § 
m nm 1 


4. Instrument Servos with Coulomb Friction and 
Stic talen 
a. A second-order servo with coulomb friction, 
(see Fig. 3-7(a)). 


The error equation for a step-displacement input is 
JE + f—F + KE +C Sgn (E) = QO 


from which the isocline equation becomes 


Pati KE e Sgn (E) 


MJ+£ Md+f 





[ze 


esi 
Il 
Ke 


Note: Let 


ne Me age 
“a “cE ac 


The second term introduces the effect of the nonlinear- 
ity. This effect is seen to be translation of the focal 


meint to E = + when —& is negative and to E= - = when 


E is positive. 


tHe KIO 


= 0 (or E axis) is a dividing line. (See 


Ed 


Fig. 3-7). The trajectory terminates at D because the 
drive torque produced by the error at D is less than the 
coulomb torque C. 

b. If stiction is present the isocline equation 
becomes 


_ _ _KE S_Sgn(E) 


MJ+£  MJ+f 





esk 


This.moves the focal point to - 2. (See Fig. 3-7(b)). 
A small initial disturbance (such as might be represented 
at point A) would not produce sufficient torque to pull the 
system free from the stiction. For a large disturbance 
(such as A') the system does pull free and begins to 
acquire a velocity. 
c. When a ramp input is applied and coulomb fric- 


tion is present but stiction may be neglected, the isocline 


equation becomes 


fw. + Cc Sgn (0, - E) 





. KE 
2 ee MI+f£ 
where W. = e. is the input ramp velocity and the coulomb 


friction is designated by C sgn (6, -) = C Sgn (6) because 
the coulomb friction effect is controlled by the relative 
velocity between the friction surfaces at the output. 

The focus is plit by the C sgn (6, -i) term, and is 
translated by both C and fw. . The foci thus appear at 
E = (fw, + C)/K; they are displaced + C/K from the point 


£wW. 
ae 


K eo 
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Trajectory 0 | 0.3 


(b) 


Fig. 3-7. Phase trajectory for an instrument 
servo with Coulomb friction and 


stiction: Step response. 
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Pie scviGdingelahe 15 a horizontal straighteline at 


E = w, 
NOEes 
6 - B= & 
ak 
=O. Ps 0, = a, 


(See Fig. 3-8(a)). For a ramp input suddenly applied at 
sO them imn tld le condltrons ane E(0)=0,E0)= w,, which 
correspond to point P on the figure. 

d. When stiction is present but discontinuous, 
the isocline system of Fig. 3-8(a) is not altered. But the 
phase trajectory cannot break free from the dividing line 

& 


at E = KR It must remain on the dividing line until E = 


(See Fig. 3-8(b)). 


An 


Crs APPLICATIONS OF THE ISOCLINE METHOD 

Since there are many methods of solution of differen- 
tial equations which are comparatively simple on the phase- 
plane but difficult or impossible otherwise, it is rather 
obvious that the phase~plane provides a very convenient 
means of analysis of extremely nonlinear first and second- 
order differential equations. It seems that the isocline 
method is a general and useful method: with which to approach 
problems by using the phase-plane method. Application of 


the isocline method to specific problems illustrates the 


usefulness of this method. 
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Fig. 3-8a. Phase portrait of an 
instrument servo with 
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Fig. 3-Sb. EFfeck fof SELCE1ONn, on 
the ramp response of a 
servo with Coulomb 


Frictven. 
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1. First-Order System 


The first-order equation has the form 


ae en) (3-1) 


Assume that the function f(x,t) is continuous and single- 
valued with the possible exception of certain singular 
points. The isocline method requires a graphical construc- 
tion performed on axes of x and t and yields a solution as 

a curve. For any given point on the x-t plane, the numeric- 
al value of £(sgpe0)e andwthus of a can be calculated. 


Example 1: Find a solution for the equation 


an 7 X 1 \e 
Consider the initial conditions x da O at ty = 0, ang 
also x, = O.2 at tO = 0. For this linear equation, the 


isoclines are straight lines given by the relation m = xt. 
A family of isoclines, each carrying its line segments of 
proper slope is shown in Fig. 3-9. Solution curves are 
sketched in the figure, starting at each of the specified 
Hiieal CoOnmeata onc. 


Example 2: RL circuit with nonlinear L 


Gia 
4b Peek) 
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Since the inductor is nonlinear, and the coefficient L 
therefore variable, it is more convenient to express the 


relationship as 


The equation of the isoclines is obtained by letting 


ae = M, some constant, and noting that Eu(t) = E for t >0O, 
then 
ya 
R 


Knowing the magnetization curve of the inductor, one 
can obtain the transient current variation of the nonlinear 
RL circuit, ivs t. (See Fig. 3-10) 
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(Example 1) The equation for a damped pendulum is: 
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Fig. 3-9. (Example 
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Eig. 3-LO0. Example (ITI-C-1-B). 
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Thus 


2 
a K oH + Sin x 
dt 
= ov 4 Kveee Soigeesc ss — 
ax 
or 
Cv. SEs ne) 
ax V 
Let 
= al 
oh, ees LO ocala 3), 
dx V 
dv 


se is the slope at any point x, v of a possible phase- 


plane trajectory. 


Now assume a number of values for ey 
av 

AS EGUdETOMN: OF yy VS me 
O ve — Sine x 

1 Ve=- + sin x 
-1 O = = sin xX 

Le 

Z V=- 3 sin xX 
-2 7 =  Sibnee 

3 i = sin x 
-—3 y= = sin x 


(See Fig. 3-11). 
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Isocline solution. 


3 ie. 
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Example 2: Nonlinear system x + xx + x = 0 


Draw a number of straight lines ¥ = Mx through the 
Origin (corresponding to different values of M) on the 
x-x plane (as in Fig. 3-12(a)). 

Plot X vs x from equation x + xx + x = 0 for con- 
stant (different) values of x (see Fig. 3-12(a)). 

With any particular straight line in view, say M = Mi, 
determine all possible points of intersections of above 
plots. The procedure is repeated for M = M,, M., Cte 
until a fair number of isocline curves are obtained on the 
phase plane. The respective slope lines are marked on 
these isoclines. Joining these slope lines gives the dif- 
ferent trajectories or the entire phase portrait (as in 
fig. 3-12(b)). 

Example 3: Consider the case of a saturating elec- 
tronic amplifier driving a motor in an instrument servo. 
(See Fig. 3-13(a)) 

The equation in the linear region is 
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Fig. 3-12(a) Example x + xx + x = 0 
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Fig. 3-12(b). Example xX + xx + x = 0 
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where E, is the value of |E | which causes saturation and MN 


is the gear ratio. 


For saturation operation, the equations are; 


ee i Cae 
ether) ee E < -E, 
rai 
R+tpe-c E 2 +E 
Gf 1 
mm 
where 
K K K 
c= = E 
NT IL 
m 


On the E vs E phase plane the isoclines are all 
straight lines; in the linear region the isoclines are 
radial straight lines emanating from a focus at the origin. 
While in the saturation region the isoclines are horizon- 
tal, the phase trajectory is constructed in the usual 
manner (as in Fags 93-13 (b))- 


Examole 4. Oscillator Wihth belay "contuser. 
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then x + x 
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Let 


x = 
then 

ys = 
where 

i 
thus 


Fig. 3-14(a). 
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iH 

Fh 

S 
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K the phase trajectory is as in 


-K Sgn(x) 


(See Fig. 3-14(b)) 


If x is positive M= - 


if x 1s negative M= - 


If u= 
uw =e & + K Sogn(x 
M4 
x = 0 or y-axis is a dividing line. 
So. auc 
y 
> Megs 
y 


If u = -K Sgn(y) 


eee x + K Sgn(y) 


_ 


y= 0 or x-axis is a dividing line. 


(See Fig. 3-14(c)) 
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EE y 1S positive M=-3 - A 
: Se ee KK 

in oy Ls Negative M= - = 
c. If u = -K Sgn(x+ty) 


Moe x -—- K Sgn(xty) 
y 


x + y = 0 1s a dividing line. 
(See Fig. 3-14(d)) 


e. If u = -K Sgn(x-y) 


M = DS aie 156 Sgn (x-y) 
sf 


X - y = 0 is a dividing line. 
(See Fig. 3-14(e)) 


Due to the changing dividing line the system became 









unstable. 
Example 5: 
u 1 
~——_$__—$fi——- a 
1 
(s+2) (st >) 
Given: 


0 ’ r=" O75 


-Sgn(x), | x | 2 Ose> 
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Since 


x(s) _ 1 = a 
Se) (s+2) (st 5) s* + 2s ar el 


Let 
x=y 
then 
wy CyYme dx | ae 
= Bs ia Sle ERS 
thus 


5 
=> ° ae 
Y 


At u= 0, the equation of isoclines is: 


At u 


- Sgn (x) 

If x is positive, the equation of isoclines is 
My = -2 Pra 

If x is negative, the equation of isoclines is 


My = - 2y = a a 


(See Fig. 3-15). 
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3. High-Order Systems 
The phase-plane, as normally used, is restricted 

to problems which can be described by first and second- 
order differential equations. This is a fairly serious re- 
striction, since many practical systems require higher-order 
differential equations for an adequate description. In the 
discussion of numerical methods for solving differential 
equations, first-order equations are studied in the begin- 
ning. Then equations of higher order are solved by replac- 
ing each Ee ea of nth order by n equations of first-order 
and solving these simultaneously. In like manner, the dis- 
cussion of graphical methods has begun by considering the 
1socline method of solving a first-order equation first, 
then second-order equations. But further consideration may 
not be possible since there is not enough information avail- 
able initially to construct isocline curves for the several 
first-order equations that would arise from a single high- 
order equation. So, if we can factor high-order equations 
into several second-order equations, these can be solved by 
phase-plane method. The difficulty is in converting a 
high-order equation into several second-order equations. 

No convenient long hand method is available, but use 
of the digital computer permits study of such phase space 


problems. 
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IV. GENERATING THE PROGRAMS TO DRAW ISOCLINES AND 
TRAJECTORIES OF SECOND-ORDER DIFFERENTIAL 


EQUATIONS 


A. xX + F(x) x + G(x) = 0 type equations (See Appendix, 


Program 1) (Where F(x) £ G(x) can be any function of x) 


SEE (x POR + G(x) = 0 
Let 
y=x 
~. Gy. dy, a& 
* “dt dx” dt 
Let 
clive 
oe M 
then 
x = My 
thus 
x + F(x) x + G(x) 
= My + F(x) y + G(x) = 0 
y=- EIT) equation of isoclines. 
Trajectories: 
ZD¢@T(1) = Z(2) 
ZD¢@T(2) = -F(Z(1))* 2(2) - G(Z(1)) 


(By using subroutine program RKLDEQ.) 
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Example 1: xX + xx + x = 0 
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we + cued Sar EQUATIONS (SEE PROGRAM 2) 


(WHERE G(x) IS ANY FUNCTION OF x) 


22 


x + x” + G(x) = 0 
Let 
x=y, 
then 
2 . Gy gx _ 
* = ax at MY 
thus | 
2 22 
xX + x” + G(x) 
2 
= My + y + G(x) = 0 
z 
y + My+ G(x) = 0 
v4 ~ 
_ <M +/M - 4G (x) ~~ 
Vai 2) ~ 
and ™ equations of 
a isoclines 
Bs a“ 
_ oM - /M"- 4G(x) - 7 
y * 2 
Trajectories: 


ZD¢T(1) = Z2(2) 


ZD¢T(2) =-Z(2) ** 2 - G(Z(1)). 
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Example 16: 
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VY-Scale = 5.00E-0l1 Units Inch. 


X-Scale 


Ge 


Big. 


Jee) 


TYPE OF EQUATIONS 


Il 
© 


x+tx|lxl+x 
(SEE PROGRAM 3) 


x =x lx l4 x =0 


In Section A and Section B the equations of isoclines 
are solved by y in terms of x. In this section we will 


solve x in terms of y to prevent imaginary values. 
‘ 


Let 
y=x, 
then 
Ehus 
lx |xX+x+x5 ly lMy ty +x =0 
e= -|y |iy ayy 
es re ly |m + 1) equation of isocline. 
Trajectories: 
ZDéT(1) = 2(2) 
ZDOT (2) SAAT) ENE 


DABS (Z(2)) 
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Vie ese || 52 eee es 6 


x+xfxl+x = My + x ly l+ x= 0 
equation of isoclines. 


Trajectories: 


ZDOT (1) ZZ) 
ZD¢T(2) =-Z(1)* DABS(Z(2)) - Z(1) 


3. £+¢xix/4+4x=0 


~ 
+ 
~ 
me 
+ 
~ 
lI 


My + y ly | + x = 0 


x = My - yly |= - y(mM + le) equation of 


isoclines. 


Trajectories: 


ZD¢T(1) = Z2(2) 


ZDéT(2) =-Z(2)* DABS(Z(2)) - Z(1). 


desis. 


= 0 


x+xIilx]+x 


Example 17: 
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Example 19: lz |X +x +x =0 


oe 
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Fig. 4-19. xX-Scale 
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D. USE OF PROGRAMS 


1. Table l. 









Type of Equation How to Set the Subroutine and 


Function Programs 



















x + F(x)x + G(x) 
(Refer to Section A) 


(1) Need to compare standard form 
to give real function F(x) and 
G(x) in each function subpro- 


gram. 


(2) Rewrite the given -equation 
x + F(x)x + G(x) = 0 into a 
system of equations in sub- 


routine program. 


Z2D@T(1) = 2(2) 
ZD¢OT(2) =-F(x)*Z(2)- G(x) 


(F(x), G(x) should change to 


z form) 













Xx + x (1) Need to give real function 





+ G(x) 


(Refer to Section B) G(x) in function subprogram. 


Same as above (2). 






lx | x +x = 0 (1) Isoclines: 





0 


X+x|lxiltex 









24 |x | + x = 0 


(Refer to Section C) 


solve x interms of y and give 






this equation in main program. 


Same as above (2). 


Table.l. 
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b. Input deck, See Table 2. 


Sequence of Cards Punch in Column 


lst Card 


2nd Card 


State ME Ores. 


4th Card 


No. of cards 
depends on 
number of dif- 
ferent ranges 

| of slope cards. 


(N+1) th Card 


(N+2) th Card 


K 
No. of cards 
depends on 
number of tra- 
jectories. 


(K+1) th Card 


Number of Cards 
of different 
ranges of slope 










‘Number of 
trajectories 


Initial time 
Final time 


Tnrtial x 


fo Cradle y 


| 


i 


| Blank 


Table 2. 
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Different ranges 


Time increment 


One significant 
number 


Column 1-10 


{Column 11-20 








Column 1-10 Ini- 
tial 
value 





11-20 
Increment 


Column 21-30 Final 
value | 


Column 














Column 


Column 


Column 
Column 
Column 


Column 


Blank 


FOR SHIFTING THE CENTER OF COORDINATES. SEE 


NODE : 


AND COMPARE PROGRAM 10. 


= 0 


cos x + sin x 


3 


Example 20). 
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Fig. 4-20. 
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E. DISCUSSION 

In ChapterII it was shown that the isocline method is 
a very useful general method of phase-plane analysis, the 
only difficulty being the labor required to get the iso- 
clines. This difficulty has been reduced by developing 
computer programs to draw both isoclines and trajectories. 

In Section A, a program works for x + F(x)x + G(x)= 0 
type equations. This is the general case, and also can be 
done by hand labor. In Section B, a program works for 
x + a + G(x) = 0 type equations. Working this kind of 
equations by hand labor, some supplementary methods are 
needed, such as Deekshatulu's or Murthy's, but use of the 
computer saves time and labor. 

In Section C, a program works for some special prob- 
lems. In this program the equation of isoclines is first 
solved for x in terms of y (Note: In Section A and B, iso- 
cline equations are solved for y in terms of x). After 
calculating the required points the plotting routine treats 
yas a function of x. The reason for this is to avoid the 
problem of complex values in the computation of y. This 
procedure is not always successful; there are at least two 
problems for which complex values apparently cannot be 
avoided. They are: 


1oeGe ss |x 


<o 
+ 
x 
x 
it 
eB. 


2. £4+x | 


Gi 


2 


For example: Equation xX + x \ x ee oe 0 
bet 

ey, 
then 

x = My 
thus 

X+x]fxil4 owe My ty ly |+ moe 
Case a: Solving for y in terms of x 


when 


y is 


Case 


when 


2 


Nee x? = 0 


-M + Jw? - 4 x2 
Z. 


va = 
M< lee 
complex. 
b: Solving for x in terms of y 
2 
x" = - [y ly | + my] 
co ily ya Moa 
Ly ly es My J 1S positive x becomes complex. 


V. GENERATING PROGRAMS TO DRAW ISOCLINES 
AND TRAJECTORIES WITH DIVIDING LINES 
OF THE SECOND-ORDER DIFFERENTIAL 


EQUATIONS 


A. SATURATION 


1. Program for Vertical Dividing Lines only 


(See Program 4) 


Refer to Chapter III, Section B-l. 
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Example 1. Trajectories: X + 0.2x + x = 0 
for Inner 


regions ) Dividing lines 


x =e 
Tsoclines: 


( See Fig. 5-1) Inner linear region: 


x + 0.2% 8 x = 0 


My + O.2y + x O 


a ~ Bs 
Y aa ioe 


Saturation regions: 


x + 022x8+ 063 = 0 


Pvc. 27 270.3 =.0 
ms gue OPES: 
4+ ~M +80.2 


Trajectories: 


Inner region 


ZD¢T (1) Z (2) 


ZDOT (2) 


-O0.2* Z(2) -Z(1) 
Saturation region 


ZDOT (1) 


Z (2) 


ZDOT(2) = -O.2* Z(2) + 0.3 


Gy n & 
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Y-Scale = 2,00E-Ol Units Inch. 


166 


= 0 
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Example 2 


fe 2 
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Dividing lines 


ee a ee aa 
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} Vie | 
4 | 1 - a 
SS —_ ores Eee 
oad -F i= 
t 4 | 
' 
i iB: 
i a \ ) 
} A “it 
* ’ 
{ 
| Ie , 
| | 7 
2 
{ { 
-_- = ee ge ere + = ee ee ee a oo = ee 
' = : 
sae x. 
: ~ 2 . ’ 
’ d i ' 
| é 4 
@ 
ne ee 
ff 
ere 





Ty; 

fn amt 
Se ne Ste 
( Tr 
fing! ; 
et 
ay Pr 
a eae 
7 el ae 
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W 

“ 

— -» 
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2. Program for Vertical Dividing Lines which can be 
Rotated. (See Program 5) 


Equation of dividing lines: 
Ax + By = +C 

Isoclines in linear region: 
g + F(x)x + G(x) = 0 


assume F(x) =1, G(x) = x 


Let y=x , then x = My 








Trajectories: 


Inner ZD¢T(1) = 2(2) 

Peguiens ZDGE(2) = -1 x Z(2) - zl) 
Saturation | 2D@PML) = 2(2) 

region Sn), = RIG Le 


liege. 


Example 3: Dividing Tames?” eat OlCy ge) Can. 


Trajeckories i eaIon2x + x i—80 











ispkeg. 1 — ore X-Scale = 2.00E-01 Units Inch, 


Y-Scale 


Z2200E—-Ol Units iIneh: 
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x - O.8y = + 0.3 


ding “anes: 
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x. 
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.ample 4 


FE 
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fé + 0.2x 
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Trajectories 
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saa] ++ OF 
(y = 0) 
M+ Oe 


Dividing lines 


Foam bes a 


= 0 


aa.4 


Trajectories 
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B. DEAD ZONE. (SEE PROGRAM 6) 
(Refer to Chapter III 


Section B-2,) 


Example 6: (See Fig. 5-6) 


Inner region: 


ITsoclines: 
%+O0.2x + 0.3 = 0 
My + Ogee, + 0.3 — 0 
M = -0.3 
M+ 0.2 


Trajectories: 
{zpeT(1) = 2(2) 


| epaeneey =-0.2% Z(2) 


External region: 
Tsoclines: 


Right and left 


th Cees oc Oe 90 


Trajectories: 


ZDOT (1) Z(2) 


ZD¢T (2) 


= 0.2 * Z(2)e— Zit) Ocee 


1h 
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Y-Scale 
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Ce IDEAL RELAY. (SEE PROGRAM 7, 8, 8A) 


For x + 0.2x + u=0 


u = - K Sgn(x) 
Switching line is x = 0 
or y axis. 
Isoclines: 


At +x XxX + 0.2x - K 


i 
© 


At -x xX + O.2x% + K 


i 
O 


Trajectories: 


At +e + 0.2% - K =O 


Z (2) 


ai 


ZD¢T (2) sue? * Z(2) +k 


At -x %¥+0.2x +K=0 


( ZD¢T (1) Z (2) 


II 


| gper(2) = - 0.2 * Z(2) - K 
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Example 72). x + Of2-cn pane 


Switching line: x 


Tsoclines: 


Right of y-axis: 


Left of y-axis: 


Dirajvec touires: 


Right ZDdT (1) 
ZDéT (2) 

Left (ZDor(1) 
\ zper (2) 


(See Fig. 5-7) 


| re 


O, 


¥. 


Il 


u = 0.3 Sen x 
Or y-axis. 

- pees 

~ M + 0.2 

= Ors 

~ M + 0.2 
Ties 

- 0.2 * Z42) + 0.3 
Z(2) 


- 0.2 * Z(2) - 0.3 


Jaa eis 


Byes fi 


X-Scale = 
Y-Scale = 
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u = - 0.3 Samco 
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Example 8(b): X + 0.02x + u = 0 
u = - 0.3 Sgn(8x-y) 
(x = y) 






ee ee et ee ee 


ee ee ee 






ee we te we re et ree ee tee ee ee 





we cr g won ‘ ~~ “ 
ee - 4 me te ee + re ee alta a —— mame = Scgtacantecn + ee oe —_ 


ee tee ay ee a eS ee ae oe eae) maa eae a 4 coerce er = x . We “a 
rie f 





-_— Se tae Pea ar oe ae Coe po a ome SS om \ < 5 x 
L pit ; : pS einer ect> th cae na Sle, Ia oe aw Ce 
he i ra ra : X, \ 
: ve r x ‘S 
/ P so ‘e . ‘\ 
‘ f ej \ ‘ * 
Thi rd ra aes xs 7 
tJ ete ~ ae oy renee fmm an nee me me eine — — ~ ht enteeen Mines ane SCs Oe ee es dy — 2 | er eye ore ——_ etme oe ee lly ere mt : aan 
¢ f ¢ \ 
. i 7 é } i y 


ew ew ew re ew —— 


oo oF 
‘ 
ey 
ee eee 
See eedientiinemikanl ep 
ay a aarp: ay 





. 
ee ee ee een te 


— ome oe gh - ewe eve we 
rd 





Pn 
we 
* 
Pe xk 
Lat oid 
fo a 
7. 
y 
4 
athat 
Pr at 
Pah 


wR re oe 


2.Q00E-Ol Units Inch. 
2.O00E-O1l Units Inch. 





dite 


eal 


- (a -_— —— 


y) 





. an a wee nee ae eee ‘= 4 = 
! 7 oo , | | : 
3 ! | | ! 
' ' ; y i | : 
| ; j | | \ : | 
: ! | : | | | | 
' ; ‘ ! | j 
; , | i 
} : { i | 
t ; i | rt 
1 2 (aa ee eee ae Stee _ 
' E ‘ . Ve 
| : : Powe | : 
: ' | : | \ ' 
' i | ' | 
t ; j | j ’ 
t H ' { ! 
. 1 ; 





4 
| 
| 
| 


2?.00E-0O1 Units Inch. 
= >) OOEmOd Unione. 


dS fa, 


t 9.02x + u = 0 


- 


-0.3 Sdn (8xty) 


Y 


om 





X-Scale 
Y-Scale 


X 


Example 8(c): 





$ 
a 
, 
. 
= { 
=e 
| 
4 
= 
I 
| 
a 


| 


- i 
r 
' 
eet ce ens ee ee 
| 
| 
U 
4 
4 
| 
1 
} 
{ 
Bee 


5-8(c). 


gaol 


Example 5-8(d): x + Q.02x + u =0 
u = - 0.3 Sgn(4x-y), x =y 
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Example 5-8(e): + 0.025 ree < 
u = - 0.3 Sgn (4xty), xX = Y 
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Example 5-8(f): xX + 0.02x + u = 0 
u = - 0.3 Sgn(2x-y), x = y 
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Di RELAY WITH DEAD ZONE (SEE PROGRAM 9) 


For ¥ +x+x+t+u=0 


U -K Sgn(x) 


Lsoclines: 
At dead zone region: 
%S+xextx=0 


My + y +x = 0 


=a 





¥"“Mdti1 
Pema eeiit = 
%+x+x -K = 0 
My + y a < -K = 90 
— Sark 
Yue id 


At lefts 


My' + yor x + K =0 
ee ae 
Y¥"“MFi1 


Trajectories: 


At dead zone region: 


ZDO6T(1) = 2(2) 
emcees ae a (O) 
Zwol (Zee) - 21) 
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Atria Gite 


; ZDéT(1) = Z2(2) 
ee ewe Geel ly 
ZDédT(2) =-2(2) - Z(1)+K 
Pea ters 
; ZO) Lie) 
met x + x + KK ae On 
ZDOE(2) == 2Z(2)- Zia 
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Example 9: 


x + 0.2 +o = 


u = - 0.3 Sgn(x) 


dividing lines: x 


2.00E-01 
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0) 
I=? 
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i il 
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Units Inch. 


2.00E-0Ol Units Inch. 





Example 10: So tases ty ea) 
U = SSCs Same) 


dividing lines: x = + 0.2 


ee 








Fig, Gr X-Scale = 2.00E-01 Units Inch. 
= 2.00E-Ol Units Inch. 


IV. CONCLUSTONS 


Chapter III summarizes eleven useful methods. Some of 
them draw direction fields to get trajectories. Some of 
them are good supplements to the isocline method. Some of 
them, by using geometrical relationships, draw the trajec- 
tories step-by-step. In general, none of these methods is 
a perfect one. For specific problems one may be better than 
the others in terms of the labor required. But if we want 
to look at the overall picture of some nonlinear system 
characteristics such as stability limit cycles, etc., none 
of them can do this easily. All need lots of time and hard 
labor. 

The isocline method has many advantages. Programs 
have been developed to compute and draw both isoclines and 
trajectories. All the programs in Chapter IV and V, first 
draw isoclines, then draw trajectories from given initial 
conditions. In analysis and design work, isoclines may not 
be needed for all problems. Conversely, sometimes only 
isoclines are needed and we do not need trajectories. Use 
of isoclines with superposed trajectories gives a global 
pitcure. That is, it indicates what can happen to trajec- 
tories for all initial conditions, and tells whether the 


system is stable or not, has limit cycle or not, 272. ete 


oe 


The programs generated in Chapter V add dividing lines 
to the isoclines and trajectories. By changing or shifting 
dividing line position, the trajectories are shaped. This 
is helpful for analysis and design work. 

In this thesis, all graphs are plotted by machine and 
all are two dimensional. The variable symbols y vs x or 
x VS9xX are used. There are no time variableswin the plot. 
If the system equation involves time as a variable, we may 
choose the 4 method to work it or apply some transforma- 
tion, changing the equation into the standard type forms. 
Tn general, machine time required to compute the isoclines 
and trajectories of a second-order nonlinear differential 
equation was about 60 sec. This isS an appreciable saving 


compared with hand labor. 
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